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ABSTRACT 

We extract Transit Timing Variation (TTV) signals for pairs of transiting planet candidates that 
are near first-order Mean Motion Resonances (MMR), using publicly available Kepler lightcurves 
(Q0-Q9). We are able to solidly confirm the planetary nature of 12 new pairs because their TTV 
signals exhibit sinusoidal variations at the expected timescale and with reasonable phases. Applying 
the analytical TTV expressions to these systems, we determine planet masses for three pairs of them 
(KOI 148, 500, 1589), based on the fact that they have nearly zero TTV phase shifts and therefore 
likely zero free eccentricities. For the remaining nine systems (KOI 152, 248, 829, 869, 877, 880, 898, 
1270, and 1336) that show non-zero TTV phase shifts, we are able to obtain their mass upper limits, 
mass ratios and thus density ratios. Combining another 11 previously confirmed Kepler near-MMR 
pairs, we find that the whole sample is generally consistent with the fiducial mass-radius relation, 
m oc i? 2 - 06 , and most of them are like twins with similar sizes and densities, although a few distinct 
exceptions, such as KOI 148, 869 and Kepler 36. In addition, there seems to be a marginal pattern, 
i.e., for slightly more pairs, the inner one has a smaller radius but larger density than that of the outer 
one. Finally, we discuss the uncertainties and implications of our findings. 
Subject headings: planetary systems— planets and satellites: detection and dynamical evolution 



1. INTRODUCTION 

As of today, Kepl er mission has found more than 2300 
planet candidates dBatalha et al.l 120121 : lOfir fc Dreizferl 
l2lMllTuan"g et al.ll2012l T and more than 880 of them are 
multiple transiting candidates (IFabrvckv et al.l l2012bl ). 
One particular interest in these multiple transiting can- 
didates is Transit Timing Variation (TTV), i.e., the de- 
viation of transit time from strict periodicity, which 
is caused by the gravitational interaction with other 
planets. By comparing the TTV modelling to the 
TTV data, one can confirm the existence of plan- 
ets and/or even constrain their masses and orbit s 
properties (jHolman fc Murray! 120051: lAgol et al.l I2005D . 
Recently, thanks to this TTV technique, many Ke- 
pler plane t systems have bee n confirmed and char- 
acterized (iHolman et al.1 120101: Lissauer et al.1 1201 lbt 

2012at ISteffen et al.1 



Cochran et al.l 120111: Ford et al 



2012t IFabrvckv et al.ll2012aHNesvornv et al.ll2012ft . Most 



of these planet systems are near Mean Motion Resonance 
(MMR). This is expected as planets in MMR build up 
their mutual interactions and thus induce particular large 
TTV signals. 

One of the keys to characterize planets with TTV is 
to find a good TTV model that can fit the TTV data. 
To do this, all the above studies rely on N-body simula - 
tions to calculate the modelling TTV (jVeras et al.ll2011l ). 
However, due to the considerable number of free param- 
eters, such fits is computationally expensive, and thus 
it is not practical for characterizing the tons of Kepler 
candidates. Moreover, the results given by N-body sim- 
ulations do not have a clear dynamical i nterpretation of 
the system. On the other side, recently, Lithwick et al. 



(|2012f) (Paper I hereafter) derived a new formulae for 
the TTV from two near-resonant planets, and proposed 
a new method to analyze the TTV data. Compared to 
the method of N-body simulation, the new method is 
much less time-consuming (hence suitable for analysis of 
a large sample), and it is straightforward to obtain the 
mass and eccentricity information of the system. 

In this paper, first, we measure and present the TTV 
of 12 new KOI pairs near first order of MMR, which 
show strong evidences to confirm their planetary nature. 
Then, applying the new method (Paper I), we analyze 
the TTV of each system, which allows us to further ex- 
tract their masses and orbital eccentricities. This pa- 
per is organized as the follows. In §2, we describe our 
method of TTV measurement and analysis. In §3, we 
present our results: the overall characteristics of the 12 
pairs and their individual properties. In §4, we discuss 
and summarize the paper. 

2. METHOD 

2.1. TTV Measurement 

The data used in this paper are the long cadence (LC), 
"corrected" light curves (PDC) of KOIs from Q0 to Q9, 
which are available at Multimission Archive at STScI 

(mastQ)- 

To measure the mid-transit times of each candidate, 
we do several irritations as described below. 

• Step 1. Segments. We first, make a guess of the 
mid-transiting times using the linear ephemcrides 
base d on the epoch and pe riod reported in the cata- 
log of lBatalha et al.l (|2012f) . Surrounding each mid- 
transit time, we cut off a segment of light curve 
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with length of 4 times transit duration. Segments 
which are overlapped by more than one transit du- 
ration with segments of other planet candidates are 
removed. Then, all the selected segments are de- 
trended using a linear based line, and their flux are 
normalized around unity. 

• Step 2. Template. Segments from last step are 
then superimposed and fitted as a whole to form 
a template using the t heoretical transiting model 
((Mandel k AgoA 12002ft . During the fitting, the 
epoch, planet-star radius ratio, transit duration, 
impact parameter, and linear limb-darkening coef- 
ficient are allowed to vary. 

• Step 3. Transit Times. Segments are fitted indi- 
vidually to the template. In each fitting, only the 
mid-transit time is allowed to vary and all other 
parameters are fixed on the values given by the 
last template. The transit time is determined by 
Levenberg-Marquardt minimization of \ 2 , and its 
unce rtainty is estimated from the covar iance ma- 
trix (|Press et al J (19921: lMarkwardtl[2009l ). The up- 
dated transit times are then used as input in step 1 
for the next irritation, and they converge generally 
in two to three irritations. 



We tested our TTV measurem ents by comparm; 
to those already pub l ished ones (|Holman et al.l 1201 



Lissauer et~al"1 l2011bt ; iCochran et all 12011 



2012a. b; Steffc n et al.ll2012t iFabrvckv et al. I l2012aft . and 



ng 
L0; 



Ford et all 



found they are generally consistent with each other. 

2.2. TTV Analysis 

In this paper, we focus on planet candidate pairs which 
are near first order MMR, i.e, 



P 



where P and P are the orbital periods of the inner and 
outer candidates, respectively. Through out this paper, 
we adopt such a convention: properties of the outer (in- 
ner) candidates arc denoted with (without) a superscript 
" ' " . The proximity to resonance is defined (see also in 
Paper I) by 



A 



P 3-1 
P 



1. 



J 



(2) 



Our method of analyzing TTV near resonance has been 
justified and presented in Paper I. Here, we briefly sum- 
marize it and describe the specific procedure used in this 
paper. 

2.2.1. TTV fitting 

As derived in Paper I, the transit time series of two 
planets near first order MMR are two linear+sinusoidal 
curves in the following forms. 



t = T + P * n + \V\sm(X j + ZV) 
t =T +P * n + \V' |sin(A J ' + ZV) 



(3) 



where, t and t , n and n , T and T , P and P and V and 
V arc the transit time, transit sequence number, transit 



offset, transit period, and TTV complex (see Eqn.7) of 
the inner and outer planets, respectively. And 

x j =j\ -U-i)x 

-(t-T) + —(j-l)(T-T) (4) 



P'/jA' 



P 



is the longitude of conjunctions, where A = ^p-{t~~T) and 
A = (t — T ) are the mean longitudes of the inner and 
outer planet, respectively. The TTV period, which we 
call super-period, is 



P 1 



i|A|' 



(5) 



and 



Attv = \V\, 
A 



A, 



W 



(6) 



are the TTV amplitude and phase, respectively. 

To fit the TTV measurements we perform the following 
procedure. 

• Step 1. Fit the transit times with linear formula, 
which give us the first guesses of T(T ) and P(P ). 

• Step 2. Calculate X 3 using the T(T') and P(P') 
obtained from last fitting. 

• Step 3. Refit the transit times using Eqn(3]to get 
new T(T'), P(P') and V(V j ). Note, in this fitting, 
all X 3 are fixed on the values obtained at step 2. 
The best fit is determined by Levenberg-Marquardt 
minimization of x 2 , and the uncertainties of fit pa- 
rameters are estimated from the cov ariance matrix 
((Press et al.lll992l : lMarkwardtll2009ft . 

• Step 4- Calculate the relative differences between 
T(T'), P(P'), P 3 and those obtained from last 
fit. If all these differences are less than 10 -3 , then 
choose the latest fit as the final result, otherwise 
goto step 2. Generally, the fitting converge in less 
than five loops. 

The direct goal of the above TTV fitting is to obtain the 
TTV amplitude and phase, as well as their error bars 
and the goodness of the fit. 

2.2.2. TTV interpretation 

The advantage of this TTV analyzing method is that 
the TTV amplitude and phase (or TTV complex) ex- 
plicitly reveal the masses and eccentricity of the system 
(Paper I), i.e., 



V = P 



V =P 



-/ 



3 Z. 



free 



2 A 



7T j A 



-9 



fix 



2 A 



(7) 



where \i and p, are the mass ratio of the inner and outer 
planets to the star, respectively, / and g are sums of 



TTV near Resonance 



3 



Laplace coefficients, as listed in Table Af of Paper I, and 



Zf lee is the complex conjugate of Zf IC 



a linear combination of the free complex eccentricities of 
the two planets. For typical Kepler system considered 
here, /i < 1CT 4 , p! < 10~ 4 , and |A| > 1CT 2 , then the 
forced eccentricities of the two planets are relative small, 
i.e., \z forced \ ~ p/\A\ < 10~ 2 and \z' forced \ ~ p /\A\ < 
10~ 2 . Thus, taking an approximation of order of unity, 
\Zf Tee \ in Eq[7]could be considered as the true eccentricity 
of the two planets. 

In principle, p, p and Z f * oe (both real and imaginary 
parts) can be inferred by inverting Eqn. [7] In reality, 

however, Z f *. cc is highly degenerated with p and/or p 
(Paper I). Nevertheless, as argued in Paper I, one can 
use the TTV phases distribution (Eqn|H|) to statistically 
estimate the masses and eccentricity. If the phases are 
close to (fittv = and <f) ttv = 180° then the system prob- 
ably have very small free eccentricity, and the masses of 
the two planets are likely to be the nominal masses solved 
from Eqnj7]by assuming \Z* free \ = 0, i.e., 



V A 



P'g 

VA 



7TJ 



Pf 



-J 2/3 (J-1) 1/3 , 



(8) 



Otherwise, if the phase arc significantly away from 



and 



180°, then it indicates large eccentricities 



of the system, and the nominal masses are likely their 
upper limits. 

Furthermore, we can estimate the mass ratio of the 
pair from Eqn|7l i.e., 




if Z lrcc < |A| 



if otherwise. 



(9) 



which allows us to further obtain their density ratio p / p. 
Note the two mass ratios in Eqn[9]are comparable to each 
other except for the case of 2:1 MMR, where they differ 
by a factor of ~ 3. 

3. RESULTS 

3.1. Overview 

We plot the TTV data for the 12 systems in Fig flM 
Their TTV generally follow the theoretical curves with 
expected period and reasonable phases, suggesting that 
they are real planets pairs near MMR. Further TTV anal- 
yses reveal the masses and orbital eccentricities of each 
system. Table [1] and table [51 respectively, summarize the 
key physical properties and the results of TTV analy- 
ses of the 12 planetary systems confirmed in this paper. 
To visualize the whole results, we plot all the systems 
in the TTV phase diagram (Fig[5]) , mass-radius diagram 
(FigEJ) , and radius-density diagram (Fig|TJ). 

3.1.1. TTV phase diagram 



We plot in Fig JS] the TTV phase of inner one vs. that 
of the outer one for each pair. One striking feature is 
that almost no system is located in quadrants II and 
Fv0, i.e., theoretically forbidden region. This is another 
independent evidence showing that these candidates are 
real interacting planets near MMR, otherwise the phase 
distribution should not follow the theoretical prediction 
but probably be random. Furthermore, the location of 
each system in the TTV phase diagram reveals the dy- 
namical information of the system. As justified in Paper 
I, generally, systems with little phase shift, i.e., close to 
(0,71-), including KOI 148, 500 and 1589, are likely to 
have zero free eccentricities and their masses are likely 
close to the nominal masses. While other systems with 
significant phase shift, including KOI 152, 248, 829, 869, 
877, 880, 898, 1270 and 1336, are likely have significant 
free eccentricities (on the order of A or even larger), and 
the nominal masses are likely their upper limits. 

3.1.2. mass-radius diagram 

We plot the nominal mass vs. radius size for each 
planet (24 in total) in Fig|51 As can be seen in this mass- 
radius diagram, most of them have radii between 2-4 R ffi 
and density around that of the Neptune and the Earth. 
Some systems are located at extremely high density re- 
gion, including 152.01, 829.01, 869.03, 1270.01. For these 
systems, such high densities should be considered only as 
upper limits. In fact, the true densities could be much 
smaller given that (1) the true masses are much lower 
than the nominal mass because of significant free eccen- 
tricities in these system, or (2) the planetary radii might 
be underestimated given their relative large uncertain- 
ties (see M4. 1 p . If the radius increases by 25% then the 
planet density will drop by a factor of 2. In addition, 
KOI 152, 500, and 1589 are systems potentially involved 
in resonance chains. Further individual analyses (i )3.2p 
show that KOI 152.03 and 152.02 are likely affected by 
the resonance chain and the masses shown here might not 
fully reflect their true masses, and the resonance chain 
effect is less important to KOI 1589.01 and 1589.02 but 
unclear to KOI 500.01 and 500.02. 

3.1.3. radius- density diagram 

We calculate the density ratio, i.e., p j p (outer over 
inner) for all the 12 pairs plus 11 pairs of previous con- 
firmed Kepler planets, and plot them vs. R /R in Fig|7l 
Ons striking feature is that all the system as a whole 
show a trend: density decreases with radius, which is 
constant with the theory of planetary internal structure 
and generally f ollows the fiducial m ass-radius relation, 
i.e., m oc i? 2 06 (jLissauer et~aT1l2011al) . Further more, for 
most systems, they are close to the origin, suggesting 
the two planets of each pair are like twins with similar 
sizes and densities. Nevertheless, there are some sys- 
tems, such as KOI 869, 1270 and Kepler 36 , with two 
planets near MMR having very different sizes and densi- 
ties. In addition, there seems to be a pattern, although 
marginally, i.e., for slightly more pairs, the inner member 
has a smaller radius but larger density than that of the 

2 Given the relative large phase uncertainty, it is reasonable that 
some of them, such as KOI 148, 877 and 880, are near the boundary 
which separate the allowed and forbidden zones. 
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outer one. Finally, wc caution that the results shown in 
this figure subject to large uncertainties, especially for 
KOI 869, 877, 898 and 1270. 

3.2. Individual Systems 

1. KOI 148.01 and 148.02 

The KOI 148 system has three transiting objects 
148.01, 148.02 and 148.03, with orbital periods 
about 4.78, 9.67 and 42.90 d and radii about 2.1, 
3.1 and 2.4 R®, respectively. The inner two are 
near 2:1 MMR with A = 0.012 and show expected 
TTV (FigHJ. Their TTV phases are near the 
boundary between allowed and forbidden zones but 
close to (0,7r), suggesting that they arc likely to 
have zero free eccentricities and their masses are 
close to the nominal mass, to ~ 14.3 and to ~ 9.2 
Mg. The outer objects, 148.03, is not close to any 
MMR and thus should have little contamination to 
the TTV of the inner pair. 

2. KOI 152.03 and 152.02 

The KOI 152 system has three transiting objects 
152.03, 152.02 and 152.01, with orbital periods 
about 13.5, 27.4 and 52.1 d and radii about 2.6, 
2.8 and 5.8 respectively. The inner two have a 
orbital ratio of 2.037 and thus near 2:1 MMR with 
A = 0.018 and an expected TTV period of 851 
day, which is clearly shown from both their TTV 
periodograms. Their TTV phases have significant 
shift away from (0, tt), suggesting they are likely to 
have significant free eccentricities (on the order of 
A at least). Combining their TTV amplitudes, we 

estimate their mass ratio to be m /to ~ 0.48, den- 
sity ratio p I p ~ 0.40, and their mass upper limits 
m nom = 65.3 and m nom = 10.1 M e . 

Note, above results ignored the effect of object 
152.01, which has a period ratio about 1.9 and 
thus A = 0.05 with 152.02. Although such an 
outer pair is much further away from 2:1 MMR 
as compared to the inner pair, it could still has a 
significant impact given the relative large radius of 
152.01. As can be seen from FigHJ, the TTV peri- 
odogram of 152.02 has a broad power distribution 
from period 500-1000 d, which is consistent with 
both super-periods of the inner (P J = 851 d) and 
outer (P J = 526 d) pairs, and thus leads to compa- 
rable fits for the TTV data 152.02 in both cases. In 
addition, the TTV periodogram of 152.01 also show 
some weak peaks around P J = 526 d. All of these 
clues lead us to conclude that the TTV of the in- 
ner pair, i.e., 152.03 and 152.02, is contaminated by 
152.01, and the results of our TTV analysis given 
in table [1] and [5] might not fully reflect the nature 
of this system. This system is particularly interest- 
ing for future studies on modelling the TTV of a 
resonant chain, which may provide further valuable 
insight into the dynamics of 3-body resonance. 

3. KOI 248.01 and 248.02 

The KOI 248 system has four transiting objects 
248.03, 248.01, 248.02 and 248.04 with orbital pe- 
riods about 2.58, 7.20, 10.91 and 18.60 d and radii 
about 1.9, 2.7, 2.5 and 2.0 R®, respectively. The 



middle two, 248.01 and 248.02 are near 3:2 MMR 
with A = 0.01 and show expected TTV (FigHJ. 
Their TTV phases (FigJSJ) have significant shift 
away from (0,7r), suggesting they are likely to have 
significant free eccentricities (on the order of A at 
least). Combining their TTV amplitudes, we es- 
timate their mass ratio to be m /»n ~ 0.60, den- 
sity ratio p I p ~ 0.73, and their mass upper lim- 
its m nom =8.9 and m nom = 6.7 M$. The other 
two smaller objects 248.03 and 248.04 are not in- 
volved in any first or second order of MMR and 
thus should have little contamination to the TTV 
of the middle pair. 

4. KOI 500.01 and 500.02 

The KOI 500 system has five transiting objects 
500.05, 500.03. 500.04, 500.01 and 500.02 with or- 
bital periods about 0.99, 3.07, 4.65, 7.05 and 9.52 
d and radii about 1.4, 1.5, 1.6, 2.6 and 2.8 R®, re- 
spectively. The outer two, 500.01 and 500.02 are 
near 4:3 MMR with A = 0.012 and show expected 
TTV (FigHJ). Their TTV phases are close to (0, tt), 
suggesting that they are likely to have zero free ec- 
centricities and their masses are close to the nomi- 
nal mass, to ~ 6.2 and m ~ 10.5 M®. 

For the inner three smaller objects, 500.05 is not 
involved in any MMR, and 500.03 and 500.04 are 
close to 3:2 MMR. In addition, 500.04 and 500.01 
are close to 3:2 MMR, and 500.04 and 500.02 are 
close to 2:1 MMR. In total, there are 4 pairs 
near MMR in the system of KOI 500, and inter- 
estingly, all of them have nearly identical super- 
period, P J ^191-192 d. Nevertheless, as seen in 
FiglU only 500.01 and 500.02 show the expected 
TTV, and the other objects seem to have little ef- 
fects (probably because their relative smaller sizes 
and thus likely smaller masses). In any case, KOI 
500 is a particularly interesting system with 4 bod- 
ies involved in 4 MMR of the same super-period. 
Future studies are needed to further unveil such a 
peculiar system. 

5. KOI 829.01 and 829.02 

The KOI 829 system has three transiting objects 
829.02, 829.01 and 829.03 with orbital periods 
about 9.75, 18.65 and 38.56 d and radii about 1.9, 
2.9 and 3.2 R®, respectively. The outer two, 829.01 
and 829.03 are near 2:1 MMR with A = 0.034 and 
show expected TTV (Fig©. Their TTV phases 
(Figj5]) have significant shift away from (0, 7r), sug- 
gesting they are likely to have significant free ec- 
centricities (on the order of A at least). Combining 
their TTV amplitudes, wc estimate their mass ra- 
tio to be to /m ~ 1.23, density ratio p / p ~ 0.93, 
and their mass upper limits m nom = 91.1 and 
m' nom = 30.5 M®. Note, in 2:1 MMR, the true 
mass of the inner one may differ a lot from its nom- 
inal mass even if the free eccentricity is small. For 
example, if \Zf ree \ = A = 0.034, then the true mass 
of 829.01 is likely smaller than its nominal mass by 
a factor of 4 (EqnJT]), i.e., to ~ 23 Mj. The inner 
smaller object, 829.02, is not close to any MMR 
and thus should have little contamination to the 



TTV near Resonance 



■5 



TTV of the middle pair. 

6. KOI 869.03 and 869.02 

The KOI 869 system has four transiting objects 
869.04, 869.01, 869.03 and 869.02 with orbital pe- 
riods about 3.22, 7.49, 17.46 and 36.28 d and radii 
about 1.8, 2.7, 2.7 and 3.2 i?®, respectively. The 
outer two, 869.03 and 869.02 are near 2:1 MMR 
with A = 0.039 and show expected TTV (Fig®. 
Their TTV phases (FigEJ) have significant shift 
away from (0, n), suggesting they arc likely to have 
significant free eccentricities (on the order of A at 
least). Combining their TTV amplitudes, we es- 
timate their mass ratio to be m jm ~ 0.28, den- 
sity ratio p I p ~ 0.16, and their mass upper limits 
m nom = 295.0 and m nom = 21.4 M ffi . 

Note, in 2:1 MMR, the true mass of the inner 
one may differ a lot from its nominal mass even 
if the free eccentricity is small. For example, if 
\Zfree\ = A = 0.039, then the true mass of 869.03 
is likely smaller than its nominal mass by a factor 
of 4 (Eqn®, i.e., m <~ 75 M®. Such a mass may 
still be too high and thus cause an unrealistically 
high density of ~ 21 gcm~ 3 for 869.03. Two possi- 
bilities or their combination may solve this density 
problem. One could be that the free eccentricity 
of the system is larger. If \Zf ree \ > 0.1 then the 
density of 869.03 will drop below ~ 8 gem -3 . The 
other possibility could be the radii of the planet 
are underestimated (see discussion in ij4.ll) . If the 
radius increase by 25%, then the density of 869.03 
will drop by a factor of 2. In any case, the density 
ratio p I ' p ~ 0.16 should be maintained (caution 
large uncertainties shown in table®. Such a large 
contrast in density is also observed recently in Ke- 
pler 36b and 36c. 

The inner two objects, 869.04 and 869.01, are not 
close to any MMR and thus has little contamina- 
tion to the TTV of the outer pair. 

7. KOI 877.01 and 877.02 

The KOI 877 system has only two transiting ob- 
jects 877.01 and 877.02 with orbital periods about 
5.95 and 12.04 d and radii about 2.4 and 2.4 i?®, 
respectively The two planets are near 2:1 MMR 
with A = 0.011 and show expected TTV (Fig®. 
Their TTV phases (Fig® have significant shift 
away from (0, 7r), suggesting they are likely to 
have significant free eccentricities (on the order of 
A). Combining their TTV amplitudes, we estimate 
their mass ratio to be m jm ~ 0.57, density ratio 
p I ' p ~ 0.60 (caution large uncertainties shown in 
table®, and their mass upper limits m nom = 21.4 

and "Com = 4 -0 M e- 

8. KOI 880.01 and 880.02 

The KOI 880 system has four transiting objects 
880.04, 880.03, 880.01 and 880.02 with orbital pe- 
riods about 2.38, 5.90, 26.44 and 51.53 d and radii 
about 1.4, 2.3, 4.0 and 5.3 i?®, respectively. The 
larger outer two, 880.01 and 880.02 are near 2:1 
MMR with A = -0.026 and show expected TTV 



(Fig®. Their TTV phases (Fig® have significant 
shift away from (0,7r), suggesting they are likely 
to have significant free eccentricities (on the order 
of A). Combining their TTV amplitudes, we es- 
timate their mass ratio to be m jm ~ 0.82, den- 
sity ratio p j p ~ 0.34, and their mass upper limits 
'mnom = 50.9 and m nom = 17.4 M®. The other 
two smaller objects 880.04 and 880.03 are not in- 
volved in any first or second order of MMR and 
thus should have little contamination to the TTV 
of the middle pair. 

9. KOI 898.01 and 898.03 

The KOI 898 system has three transiting objects 
898.02, 898.01 and 898.03 with orbital periods 
about 5.17, 9.77 and 20.09 d and radii about 2.2, 
2.8 and 2.4 R e , respectively. The outer two, 898.01 
and 898.03 are near 2:1 MMR with A = 0.028 and 
show expected TTV (Fig®. Their TTV phases 
(Fig® have significant shift away from (0,7r), sug- 
gesting they are likely to have significant free ec- 
centricities (on the order of A at least). Combining 
their TTV amplitudes, we estimate their mass ratio 
to be m jm ~ 1.11, density ratio p / p ~ 1.92 (cau- 
tion large uncertainties shown in table®, and their 
mass upper limits m nom = 44.7 and m nom = 14.3 
M®. Note, in 2:1 MMR, the true mass of the inner 
one may differ a lot from its nominal mass even 
if the free eccentricity is small. For example, if 
\Z free \ = A = 0.028, then the true mass of 898.01 
is likely smaller than its nominal mass by a factor 
of 4 (Eqn®, i.e., m ~ 11 M®. The other inner ob- 
jects 898.02 is not close to any first or second order 
of MMR and thus should have little contamination 
to the TTV of the outer pair. 

10. KOI 1270.01 and 1270.02 

The KOI 1270 system has only two transiting 
objects 1270.01 and 1270.02 with orbital periods 
about 5.73 and 11.61 d and radii about 2.2 and 
1.6 R®, respectively. The two planets are near 2:1 
MMR with A = 0.013 and show expected TTV 
(Fig®. Their TTV phases (Fig® have significant 
shift away from (0,7r), suggesting they are likely 
to have significant free eccentricities (on the order 
of A at least). Combining their TTV amplitudes, 
we estimate their mass ratio to be m jm ~ 0.14, 

density ratio p j ' p ~ 0.39 (caution large uncertain- 
ties shown in table®, and their mass upper limits 

rrinom = 132.8 and m nom = 6.0 M®. Note, in 2:1 
MMR, the true mass of the inner one may differ a 
lot from its nominal mass even if the free eccentric- 
ity is small. For example, if \Zf ree \ = A = 0.013, 
then the true mass of 1270.01 is likely smaller than 
its nominal mass by a factor of 4 (Eqn®, i.e., 
m ~ 33 M®. Such a mass may still be too high 
and thus cause an unrealistically high density of 
~ 17 gcm~ 3 for 1270.01. Two possibilities or their 
combination may solve this density problem. One 
could be that the free eccentricity of the system is 
larger. If \Z free \ > 0.03 then the density of 1270.01 
will drop below ~ 8 gem -3 . The other possibility 
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could be the radii of the planet are underestimated 
(see discussion in £ )4.1I) . If the radius increase by 
25%, then the density of 1270.01 will drop by a fac- 
tor of 2. Indeed, we note there are quite large un- 
certainties on the planet radii of system KOI 1270 
(see table [TJ 

11. KOI 1336.01 and 1336.02 

The KOI 1336 system has only two transiting 
objects 1336.01 and 1336.02 with orbital periods 
about 10.22 and 15.57 d and radii about 2.8 and 
2.9 R®, respectively. The two planets are near 3:2 
MMR with A = 0.016 and show expected TTV 
(FiggJ). Their TTV phases (Fig© have significant 
shift away from (0, n), suggesting they are likely to 
have significant free eccentricities (on the order of 
A at least). Combining their TTV amplitudes, we 

estimate their mass ratio to be m /m ~ 0.89, den- 
sity ratio p I p <~ 0.82, and their mass upper limits 
m nom = 23.9 and m nom = 26.4 M e . 

12. KOI 1589.01 and 1589.02 

The KOI 1589 system has five transiting objects 
1589.04, 1589.01, 1589.02, 1589.03 and 1589.05 
with orbital periods about 4.22, 8.73, 12.88, 27.43 
and 44.55 d and radii about 1.3, 2.2, 2.4, 2.4, and 
2.0 R®, respectively. 1589.01 and 1589.02 are near 
3:2 MMR with A = -0.016 and show expected 
TTV (FiggJ). Their TTV phases are close to (0, tt), 
suggesting that they are likely to have zero free ec- 
centricities and their masses are close to the nomi- 
nal mass, m ~ 17.6 and m ~ 16.5 M®. 

For the other three objects, the outer two, 1589.03 
and 1589.05, are not close to any first or second 
order of MMR and should have little contamina- 
tion to the TTV of 1589.01 and 1589.02. The in- 
nermost one, 1589.04, is close to 2:1 MMR with 
1589.01. Such a pair have a A = 0.033 and super- 
period of P J = 133 d. Considering the relative 
large A (proximity to resonance, which is two times 
larger than that of the pair 1589.01 and 1589.02), 
and the relative smaller radius of 1589.04 as com- 
pared to 1589.01, we expect that 1589.04 should 
have little contamination to the TTV of 1589.01 
and 1589.02. Indeed, there is no evident signal 
in the TTV of 1589.01 with period close to the 
super-period (P J = 133 d) of the pair 1589.04 and 
1589.01. 

4. SUMMARY AND DISCUSSION 

A large number of transiting exoplanet candidates have 
been, and would continually be found by the Kepler mis- 
sion, raising the demand for an efficient method to con- 
firm and characterize the planetary nature of these can- 
didates. Measuring and analyzing TTV is a promising 
technique to fill the above need, and it has been demon- 
strated to be highly successful, espec ially for those multi- 
ple transiting candida t es near MMRdHolman et al.H2010t 
Lissauer et al.l l2011bfc [Cochran et al.l 120111: iFord et al.l 
2012al iFabrvckT et al. l2012aHSteffen et alJI2012h . How- 
ever, all these studies, rely on fitting the TTV data to di- 
rect N-body simulations, which is computationaly time- 
consuming and thus more suitable to detail analysis of 
individual systems. 



In this paper, we extract TTV signals for 12 sys- 
tems of transiting planet candidates that are near first- 
order MMR, using publicly available Kepler lightcurves 
through Q9. Their TTV (FigJUgJ) generally follow the 
theoretical curves with expected period and reasonable 
phases, suggesting that they are real planets pairs near 
MMR. Based on the new analytical formulae for the TTV 
from two near-resonant planets (Paper I) , we are able to 
develop an efficient method for analyzing the whole sam- 
ple of TTV data. The direct results (tabled]) of this TTV 
analysis progress are the TTV amplitude and phase (or 
called TTV complex together, defined in the TTV formu- 
lae at §2.3 and/or Paper I), which reveal the information 
of the mass and eccentricity of the system (see Figj5][7] 
and table [T|) . We summarize our major findings in the 
following. 

• Three systems KOI 148, 500 and 1589 are likely to 
have zero free eccentricities and their masses are 
likely close to the nominal mass (EqJS]). For the 
other nine systems with significant free eccentric- 
ities (on the order of A at least), we derive their 
mass upper limits, mass ratios and density ratio 
(caution the large error bars). 

• Three systems, KOI 152, 500, and 1589, are close to 
resonance chains. We find that pair KOI 152.03-02 
are likely affected by the resonance chain induced 
by 152.01, therefore the results of their TTV anal- 
ysis shown here might probably be contaminated. 
The resonance chain effect is found to be unlikely 
to affect the TTV analyses of KOI 1589.01-02 and 
KOI 500.01-02. 

• Combining another 11 previously confirmed near- 
MMR KOI pairs, we find (1) the whole sample is 
generally consisten t with the fiducial m ass-radius 
relation m oc i? 2 06 (|Lissauer et al.l20liaj) , (2) most 
pairs are like twins with similar sizes and densities 
(caution the large error bars, see N4.ip . although 
some pairs with dramatic differences in sizes and 
densities of their members, and (3) a marginal pat- 
tern, i.e., for slightly more pairs, the inner member 
has a smaller radius but larger density than that 
of the outer one. All these clues provide valuable 
insights into the origin of such near-MMR parirs 
(see aOl . 

4.1. Uncertainties of Stellar and Planetary Radii 

Planetary radius is determined by the product of two 
independent quantities: planet-star radius ratio, which 
is based on the modelling transit light curves, and stel- 
lar radius, which is based on our knowledge of the host 
star. The radii of the planets confirmed in this pap er are 
adopted from the catalog of iBatalha et al.l (|2012[ ). and 
their uncertainties plotted in Fig(S] and table Q] only re- 
flect the uncertainties on individual fits of transit light 
curves. We caution the true uncertainties are likely un- 
derestimated given the consideration of the following as- 
pects. 

• Systematic uncertainty of planet-star radius ratio 
could be induced by the different methods and 
models used to analyze the transit light curves, 
such as the specific way for detrending and folding 
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the light curves and the specific limb darkening law 
used in fitting the transit light curves. For a sim- 
ple diagnosis, we compared planet-star radius ratio 
measurements of 6 Kepler confirmed p lanets from 
the KOI catalog (jBatalha et al.l l2012f ) and those 
reported in papers which individually claimed the 
confirmation of such planets. As can be seen in 
the middle panel of FigfSJ the standard deviation 
of such two kinds of measurements is 12.6%, which 
may reflect the order of magnitude of the system- 
atic uncertainty of planet-star radius ratio. 

• Uncertainties of stellar radii are not given in the 
KOI catalog (jBatalha et al.ll2012D . but can be gen- 
erally estimated by comparing the stellar radii 
given but different studies. We perform such an 
analysis to 37 KOIs in the top pan el of Fig [8] As 
can b e seen, the stellar radii from (Batalha et al. 
120121 ) tend to be systematically underestimated by 
~ 7% and has a standard deviation of ~ 28% as 
compared to those given by individual confirming 
papers. 

In the bottom panel of FigEl we plot the deduced plan- 
etary radii of 60 Kepler confirmed planets from individ- 
ual confirming paper s vs. those from the KOI catalog 
(jBatalha et al.l I2012D . According to principle of error 
propagation, we expect that the uncertainty of the plan- 
etary radii should be the sum of both uncertainties of 
planet-star radius ratio and stellar radius. As can be 
seen, the planetary radius deviation of 29%, although 
lower than the above expectation by ~ 12% 0, could still 
cause a large deviation in planet density of a factor of 
~ 2. Future high-precision spectroscopic follow-up ob- 
servations which can help accurately determine stellar 
radius are required to better characterize these KOI sys- 
tems. 

4.2. Origin of the near-MMR pairs 

One striking feature of the near-MMR KOI pairs is 
the asymmetric distribution around the nominal res- 
onance center, i.e., most of them have a period ra- 
tio which is a f ew percents larger th an the exact 
commensurabilitv. (jFabrvckv et al.1 l2012b[ K Such a fea- 
ture is also seen in the 12 planets confirmed in this paper, 
among which 2 of them with negative A while the other 
10 with positive A. Recently, it has been suggested that 
such an asymmetry might be induced by some dissipation 
processes, e.g., tides, which damp the eccentricities of 
planets in MMR and lead to a slow divergence away from 
exact commensurabilitv (iTerauem fc Papaloizo'ul [20071 : 
iBatalha et aI|[20Tl iLithwick fc Wul 12012ft . One predic- 
tion of such an dissipation scenario is that the near-MMR 
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3 This might be caused by that either planet-star ratio and stellar 
radius are not fully independent of each other or the samples (37 
stars and 60 planets) studied here are not large enough 



pairs should have zero-free eccentricities because of ec- 
centricity damping. Although this seems to be confirmed 
by some systems (KOI 148, 500 and 1529), we caution 
that all the other systems have small but significantly fi- 
nite free eccentricities on a order of a few percents (sim- 
ilar results were also abstained in Paper I). This may 
imply either planets could reacquire free eccentricities af- 
ter dissipation via some certain processes or maybe the 
asymmetric feature was caused by some other unknown 
mechanisms. 

More clues might be revealed from the radius-density 
diagram (Fig [7]). As can be seen, most paris are clus- 
tered near the origin where p — p and R = R, which 
might indicate that most of them are randomly paired 
together with similar sizes an d densities. Such a re - 
sult is consistent with that of iFabrycky et ail (j2012b| ). 
which finds that all paris in all KOI systems are likely 
drawn from a smooth distribution of period ratios, and 
they do not necessarily point to differential migration. 
Altho ugh these clues seems t o support an insitu sce- 
nario ([Hansen fc Mu rray 20 fH) for the formation of those 
pairs, we caution that an non-insitu scenario associ- 
ated with migration, if not dominant, might still be 
suitable for explaining the marginal radius-density pat- 
tern, namely more pairs with the inner member smaller 
but denser than the outer one, such as Kepler 36 and 
KOI 869. The general explanation from the migration 
scenario could be that the currently near-resonant pair 
were not formed near each other, but widely separated 
in the protoplanetary disk. As the disk is cooler and 
hill radius is larger (thus larger feeding zone when the 
planet was formed) outside, the outer planet accreted 
more volatile material and finally became larger but less 
dense than the inner one. After their formation, interac- 
tions with the disk c ould then enable them to migrate to 
the present location (iLee fc Peald[2002l: iZhou et al.ll2005t 
ITerauem fc Papaloizoull2007D . 

The current results shown in this paper arc subject to 
relatively large uncertainties both on planet mass and ra- 
dius. Further TTV data and high-precision spectroscopic 
follow-up observations of the host stars in the future will 
reduce such error bars, thus allow a better characteriza- 
tion to each system and provide a clearer picture on the 
origin of the near-MMR pairs. 

JWX acknowledges support from the Ontario govern- 
ment, the National Natural Science Foundation of China 
(Nos. 10833001 and 10925313), Fundamental Research 
Funds for the Central Universities (No.1112020102), Re- 
search Fund for the Doctoral Program of Higher Educa- 
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Fig. 1. — TTV evidences of four candidate paris. For each of them, we plot the best-fit theoretical curves on top of the TTV data and the 
TTV pcriodogram. In the TTV fitting panels, The vertical dashed lines denote the times when the longitude of conjunction points at the 
observer, i.e., = 0. In the pcriodogram panel, the red vertical dashed line denotes the theorectically predicted period of the TTV, and 
the two other vertical lines denotes the periods corresponding to the Nyquist critical frequencies of the two planets, respectively. The TTV 
fitting results: amplitude (Attv) and its normalized error, phase (rpttv) and its absolute error, planet radius (R), nominal mass (m nom ), 
and orbital period (P) are also shown in each TTV fitting panel. 
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Fig. 2. — Similar to Fig[T] but for another four pairs. 
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Fig. 3. — Similar to Fig[T] but for another four pairs. 
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Similar to Fig[T] but for another four pairs. 
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each panel, we print the average ± standard deviation of the ratio between values Y and X. Such analyses provide a simple diagnosis to 
the uncertainties of stellar radius, planet-star radius ratio and planetary radius. 
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TABLE 1 

Key Properties of Planets and Stars of 12 Systems. 



koi a id id' 


P P' 


R b R' b 


c ' c 


m /rn d p / p e 


Mj R^ 


_ 


d d 


R® ft® 




- 


M R 


148 1 2 


4.778 9.674 


2.1(1±0.01) 3.1(1±0.01) 


14.3(1±0.24) 9.2(1±0.27) 


0.64(1±0.52) 0.20(1±0.60) 


0.89 0.89 


152 m + 3 2 


13.485 27.403 


2.6(1±0.01) 2.8(1±0.01) 


65.3(1±0.19) 10.1(1±0.33) 


0.48(1±0.52) 0.40(1±0.59) 


1.07 0.99 


248 1 2 


7.204 10.913 


2.7(1±0.04) 2.5(1±0.05) 


8.9(1±0.08) 6.7(1±0.16) 


0.60(1±0.24) 0.73(1±0.50) 


0.55 0.53 


500 m ~ 1 2 


7.053 9.522 


2.6(1±0.03) 2.8(1±0.03) 


6.2(1±0.22) 10.5(1±0.14) 


1.70(1±0.36) 1.44(1±0.55) 


0.72 0.65 


829 1 3 


18.649 38.557 


2.9(1±0.04) 3.2(1±0.04) 


91.1(1±0.21) 30.5(1±0.21) 


1.23(1±0.42) 0.93(1±0.64) 


0.97 0.89 


869 3 2 


17.461 36.275 


2.7(1±0.09) 3.2(1±0.09) 


295.0(1±0.20) 21.4(1±0.36) 


0.28(1±0.56) 0.16(1±1.09) 


0.79 0.82 


877 1 2 


5.955 12.040 


2.4(1±0.15) 2.4(1±0.15) 


21.4(1±0.25) 4.0(1±0.27) 


0.57(1±0.52) 0.60(1±1.41) 


0.63 0.59 


880 1 2 


26.444 51.535 


4.0(1±0.01) 5.3(1±0.01) 


50.9(1±0.36) 17.4(1±0.28) 


0.82(1±0.65) 0.34(1±0.70) 


0.93 0.91 


898 1 3 


9.770 20.090 


2.8(1±0.14) 2.4(1±0.14) 


44.7(1±0.43) 14.3(1±0.29) 


1.11(1±0.72) 1.92(1±1.55) 


0.66 0.61 


1270 1 2 


5.729 11.609 


2.2(1±0.24) 1.6(1±0.24) 


132.8(1±0.08) 6.0(1±0.36) 


0.14(1±0.44) 0.39(1±1.86) 


0.83 0.73 


1336 1 2 


10.219 15.574 


2.8(1±0.04) 2.9(1±0.04) 


23.9(1±0.28) 26.4(1±0.35) 


0.89(1±0.63) 0.82(1±0.88) 


0.95 1.03 


1589 m " 1 2 


8.726 12.883 


2.2(1±0.03) 2.4(1±0.03) 


29.2(1±0.26) 19.9(1±0.35) 


0.68(1±0.61) 0.58(1±0.77) 


0.96 0.99 



a Superscirpt "m+" denotes that TTV analysis of KOI 152.03-02 is likely contaminated by the resonance chain induced by KOI 
152.01 and the masses shown here might not fully reflect their true masses and "m-" denotes that the resonance chain effect is less 
important to the TTV analyses of KOI 500.01-02 and KOI 1589.01-02. 

6 Adopted from Batalha ct al. (2012), and the uncertainties are only the error bars of planet-star radius ratio. 

c Deduced from Eqn|8] which are likely the true mass for system KOI 148, 500 and 1589 but upper limits for others. 

d Deduced from Eqn|9] and each error bar is computed as the sum of TTV amplitude errors of the pair. Note this ratio could differ a lot 
compared to the ratio of nominal masses, especially for the case of 2:1 MMR. 

e Deduced from the above mass ratios and radii. Note density ratio is independent on the uncertainty of stellar radius. 
f Deduced from the stellar rad i i and surface gravity reported from[Batallia et al. ( 2012T) . 
f Adopted from lBatalha et al.l GWFZ ). 



TABLE 2 

Results of TTV Analyses for 12 Pairs of Planets 



KOI" 


id 


id' 


j 


P 


P 


A 


Pi 


Attv 




</>ttv 


</>ttv 


Xrcd 


/ 










d 


d 




d 


d 


d 


deg 


deg 






148 


1 


2 


2 


4.778 


9.674 


0.012 


392.0 


0.0028(1±0.27) 


0.0023(1±0.24) 


-23.5±17.6 


199.5±18.5 


1.4 


1.2 


152 m + 


3 


2 


2 


13.485 


27.403 


0.016 


851.9 


0.0055(1±0.33) 


0.0182(1±0.19) 


-29.6±28.6 


123. 2± 8.7 


1.2 


1.1 


248 


1 


2 


3 


7.204 


10.913 


0.010 


367.7 


0.0064(1±0.16) 


0.0139(1±0.08) 


54.1± 9.3 


218. 0± 5.5 


1.3 


0.9 


500 m - 


1 


2 


4 


7.053 


9.522 


0.012 


191.1 


0.0059(1±0.14) 


0.0049(1±0.22) 


-9.1± 9.0 


178.0±13.3 


1.3 


1.7 


829 


1 


3 


2 


18.649 


38.557 


0.034 


571.4 


0.0116(1±0.21) 


0.0155(1±0.21) 


-33.9±11.1 


135.9±13.0 


1.2 


1.1 


869 


3 


2 


2 


17.461 


36.275 


0.039 


468.0 


0.0081(1±0.36) 


0.0476(1±0.20) 


24.0±18.6 


236.2±10.1 


1.0 


1.9 


877 


1 


2 


2 


5.955 


12.040 


0.011 


551.2 


0.0024(1±0.27) 


0.0069(1±0.25) 


-3.3±16.5 


269.8±11.4 


1.2 


1.4 


880 


1 


2 


2 


26.444 


51.535 


-0.026 


1007.0 


0.0145(1±0.28) 


0.0275(1±0.36) 


1.5±13.5 


337.2±12.3 


1.3 


3.2 


898 


1 


3 


2 


9.770 


20.090 


0.028 


357.5 


0.0051(1±0.29) 


0.0075(1±0.43) 


63.8±15.0 


213.2±23.3 


1.2 


1.2 


1270 


1 


2 


2 


5.729 


11.609 


0.013 


441.6 


0.0022(1±0.36) 


0.0258(1±0.08) 


96.6±18.5 


281. 7± 4.7 


1.4 


1.1 


1336 


1 


2 


3 


10.219 


15.574 


0.016 


323.4 


0.0125(1±0.35) 


0.0187(1±0.28) 


-61.7±17.3 


125.6±15.8 


1.1 


1.4 


1589 m " 


1 


2 


3 


8.726 


12.883 


-0.016 


272.7 


0.0089(1±0.35) 


0.0204(1±0.26) 


21.3±26.3 


184.6±14.1 


1.2 


1.2 



a superscripts "m+" and "m-" have the same meaning as in tabic [T] 



The TTV data are available at www. 
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